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Muon-induced neutrons do not explain the DAMA data 
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We present an accurate model of the muon-induced background in the DAM A/LIBRA experi¬ 
ment. Our work challenges proposed mechanisms which seek to explain the observed DAMA signal 
modulation with muon-induced backgrounds. Muon generation and transport are performed using 
the MUSIC/MUSUN code, and subsequent interactions in the vicinity of the DAMA detector cav¬ 
ern are simulated with Geant4. We estimate the total muon-induced neutron flux in the detector 
cavern to be = 1.0 x 10“® cm“^ s“^. We predict 3.49 x 10“® counts/day/kg/keV, which accounts 
for less than 0.3% of the DAMA signal modulation amplitude. 


I. INTRODUCTION 

The DAMA/LIBRA experiment [T] is a highly radiop- 
ure Nal(Tl) scintillation detector located at the Gran 
Sasso National Laboratory (LNGS) which aims to mea¬ 
sure the annual modulation signature of dark matter par¬ 
ticles [H 13]. Both the DAMA/LIBRA experiment and 
the first generation DAMA/Nal experiment reported the 
observation of an approximately annual variation in the 
number of events observed in the 2-6 keV energy range 
with a combined significance of approximately 9.3 cr [Ij. 
If the observed signal modulation is to be explained by 
the elastic scattering of dark matter particles, it would 
require that the interaction cross-section of dark matter 
with nucleons, and the mass of dark matter particles, 
to be within values that are already excluded by other 
experiments [Hltllj. 

One mechanism that has been proposed in order to ex¬ 
plain the DAMA signal modulation is the production of 
neutrons due to the scattering of cosmogenic muons in 
the material surrounding the detector [T2|-[l4|. The cos¬ 
mogenic muon-induced-neutron flux is expected to 
have an annual variation related to the mean air tem¬ 
perature above the surface of the Earth that affects the 
muon flux <I>'^ at the surface, and hence underground. 
This proposal has been disputed for a number of rea¬ 
sons [3|, notably as the annual variation of cosmogenic 
muons is approximately 30 days out of phase with the 
DAMA signal [TSHTB] . 

An extension to this mechanism has also been pro¬ 
posed, which introduces the possibility of a contribution 
to the total neutron flux from the interactions of so¬ 
lar neutrinos m- The solar neutrino-induced neutron 
flux is also expected to have an annual modulation, 
due the eccentricity of the Earth’s orbit about the Sun. 
It is shown that the phase of <I>J^ can shift the phase of 
the total neutron flux relative to 

One should note that most explanations focus on the 
phase, rather than the amplitude, of the modulation. 
A rough estimate of the modulated rate of muon-induced 
neutrons R!f^ in the DAMA/LIBRA experiment can be 
calculated as: 

~ 4.6 X 10-® events / day / kg (I) 
m 

where is taken from previous estimates at LNGS 


[ini ini, is the amplitude of the muon flux modu¬ 
lation [311 [13| (relative to $^), A is the active surface 
area perpendicular to of one DAMA/LIBRA detec¬ 
tor module, t is the number of seconds in a day and m is 
the active mass of one detector module. This rate is three 
orders of magnitude less than the modulated DAMA sig¬ 
nal, even before the acceptance for muon-induced events 
in the DAMA analysis is considered. 

Despite the latter calculation, and previous work [31 
HZlIISj, the recent paper m claims that the DAMA sig¬ 
nal can partly be explained by muons. The estimate 
presented by Davis in Ref. m. includes only a very 
approximate calculation of the amplitude of muon- and 
neutrino-induced signals. The calculation of contra¬ 
dicts the estimates from Refs. [31III1IIH1- The calcula¬ 
tion neglects that most neutrons will be accompanied a 
showering muon, and would therefore not be accepted by 
DAMA/LIBRA. Also, the value of the mean free path 
for neutrons in rock that is used (taken from Ref. [24 ] ). 
is actually for liquid argon; which has approximately half 
of the density of LNGS rock. Ref. [IS] has additionally 
shown that the required 4)5( is at least six orders of mag¬ 
nitude too large. 

It is evident from such contradictions that a full Monte 
Garlo simulation of the muon-induced background in the 
DAMA/LIBRA experiment is required. This would be 
able to model $((, the detector response and the DAMA 
event acceptance, and also any enhancement of due 
the high-Z shielding used by DAMA/LIBRA. 

In this paper, we present an accurate calculation of $(( 
and the total muon-induced background in order to show 
that the DAMA signal modulation cannot be explained 
by any muon-induced mechanism. We perform a full 
simulation of the DAMA/LIBRA apparatus and detector 
shielding in Geant4.9.6 [IS], with muons transported to 
the DAMA cavern in LNGS using the MUSIC/MUSUN 
code [27l [28] . In total, we simulate twenty years of muon- 
induced data. 


II. SIMULATION FRAMEWORK 

We perform the simulation of particle propagation in 
two stages. In the first stage, only muon transporta¬ 
tion from the surface of the Earth down to an under- 
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ground site is considered and secondary particles [5^ are 
neglected. In the second stage, the transport and in¬ 
teractions of all particles (including secondary particles) 
are fully simulated through the material surrounding the 
DAMA/LIBRA apparatus. 

A. Muon transport simulation 

The first stage of the simulation is performed using 
the MUSIC muon transport code [271 [2H]. The MUSIC 
code propagates muons from the surface of the Earth 
through a uniform rock of density p = 2.71 g cm“^ 
and records energy distributions of muons at different 
depths. The MUSUN code [27] calculates muon spec¬ 
tra from the modified Gaisser’s parameterisation that 
takes into account the curvature of the Earth and muon 
lifetime, convoluted with the slant depth distribution at 
LNGS. This parameterisation has been previously shown 
to have a good fit to LVD data [27|. The MUSUN code 
subsequently samples muons on the surface of a cuboid 
with a height of 35 m and perpendicular dimensions of 
20 m X 40 m. This cuboid includes most of the corridor 
where the DAMA/LIBRA experiment is located and a 
few meters of rock around it. 


B. DAMA/LIBRA detector simulation 

The second stage of the simulation is performed using 
Geant4.9.6 [IS]. The Geant4.9.6 shielding physics 
list has been used, and we include the muon-nuclear inter¬ 
action process. The interactions of low-energy neutrons 
(< 20 MeV) are described by high-precision data-driven 
models [30]. Previous studies [2I][3TH31| have validated 
the simulation of neutron production, transport and de¬ 
tection against data. The level of agreement is better 
than a factor of two. 

In this phase of the simulation, all primary and sec¬ 
ondary particles are transported from the surface of the 
cuboid until all surviving particles have propagated out¬ 
side of the cuboid volume. The cuboid is modeled as 
LNGS rock with a density p = 2.71 g cm“^ and a chem¬ 
ical composition as described in Ref. [35|. A corridor 
(‘cavern’) is positioned within the cuboid, such that there 
is 10 m of LNGS rock overburden, and otherwise 5 m of 
LNGS rock surrounding the cavern walls and floor. 

The DAMA/LIBRA detector housing is placed halfway 
along the length of the cavern, adjacent to a cavern wall. 
The housing is composed of LNGS concrete with den¬ 
sity p = 2.50 g cm“^ and a chemical composition as 
described in Ref. [35]. The DAMA/LIBRA apparatus 
and detector housing are described in Ref. [T] . There are 
a number of concentric layers of shielding surrounding 
the DAMA/LIBRA detector. Extending outwards from 
the detector, we model 10 cm of copper, 15 cm of lead, 
1.5 mm of cadmium, 50 cm of polyethylene and 1 m of 
LNGS concrete. 


We model each of the 25 DAMA/LIBRA detector 
modules, containing a central cuboidal crystal composed 
of Nal, in addition to light-guides and photomultiplier 
tubes [T]. The dimensions of each module, including 
a further 2 mm of copper shielding, are 10.6 x 10.6 x 
66.2 cm^. The 25 modules are placed in a 5 x 5 arrange¬ 
ment in the vertical and width dimensions of the cavern. 


III. THE MUON-INDUCED NEUTRON FLUX 

In the stage of simulation described in Section |IIB[ 
muons and secondary particles are transported through 
10 m of LNGS rock above the DAMA cavern, and also 
through 5 m on the sides and underside of the cavern. 
Integrating over the surface area of the cavern, our sim¬ 
ulation predicts = 1.0 x 10“® cm“^ s“^, excluding 
back-scattering [33] . 

In Table |lj we compare our result to the simula¬ 
tion of Wulandari et al. [20], which is performed using 
ELUKA [37] , and of Persiani [^ , which is performed us¬ 
ing Geant4.9.3 and MUSIC/MUSUN. Integrating over 
all neutron energies, our results are consistent within 
about 30% of the previous estimates. 

We additionally demonstrate a dependance of $(( on 
the dimensions of the cavern by scaling <!>(( to the cavern 
proportions used by Wulandari et al. (compare the first 
two rows in Table |l]) . We attribute this to the different 
fluxes and energy spectra of vertical and inclined muons. 

High-Z materials in the detector shielding (lead and 
copper) will lead to an enhancement of $(( which could, 
potentially, contribute to the modulated signal. Eigure[T] 
shows as a function of neutron energy, as predicted in 
the LNGS cavern and after all particles are propagated 
through the various layers of the DAMA/LIBRA shield¬ 
ing. It is shown that $(( increases by a factor of > 5 
due to the shielding. As we will discuss in Section |IV] 
this enhancement of $(( is still insufficient to explain the 
DAMA data. 


TABLE 1. A comparison of (in units of 10“^° cm“'^ s“^) 
predicted by this study, Wulandari et al. [^ and Persiani |21| . 
The column titled ‘Cavern’ indicates the three distinct cavern 
geometries used: (1) in this study; (2) by Wulandari et al. and 
(3) by Persiani. The range of considered neutron energies is 
shown, and ‘(*)’ indicates that back-scattered neutrons are 
included. 


Cavern > 0 MeV > 1 MeV > 1 MeV (*) 


This study 

(1) 

10 

4.0 

5.0 

This study 

(2) 

7.6 

5.8 

10 

Wulandari et al. 

(2) 

No data 

4.3 

8.5 

Persiani 

(3) 

7.2 

2.7 

No data 
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FIG. 1. The distribution of a function of neutron energy 
for neutrons entering the cavern in which DAM A/LIBRA is 
situated (solid lines) and entering the DAM A/LIBRA detec¬ 
tor modules after all shielding is traversed (dashed line). The 
distributions excluding and including back-scattered neutrons 
are shown in black and red respectively. 


IV. ANALYSIS 

In our simulation, each detector module is treated in¬ 
dependently and all information due to an energy de¬ 
position from any particle in the Nal crystal volumes is 
recorded. DAMA categorises events as being single-hit 
or multiple-hit if the event has an energy deposit in only 
a single crystal or in multiple crystals respectively. The 
DAMA signal region, in which the modulated signal is 
observed, is then defined only for single-hit events with 
a total energy deposit (Euep) in the range 2-6 keV. 

In the following sections, we will show that the num¬ 
ber of muon-induced events entering the DAMA/LIBRA 
signal region is too low to explain the signal modulation. 


A. Resolution and qnenching factors 


FIG. 2. The energy spectrum for all crystals with energy de¬ 
positions, with resolution and nuclear recoil quenching factors 
applied. The equivalent of twenty years of muon-induced data 
is presented. 



FIG. 3. The distribution of the hit multiplicity in events with 
i?Dep > 2 keV per crystal. At least one crystal has a total 
energy absorption in the range 2-6 keV. The equivalent of 
twenty years of muon-induced data is presented. 


these events are single-hit events, which suppresses any 
enhancement of due to interactions in the detector 
shielding. 


We model the DAMA/LIBRA experimental resolution 
by applying a Gaussian smearing to the sum of all en¬ 
ergy deposited in each crystal, using resolution parame¬ 
ters provided by Ref. [T]. As Geant4 does not account 
for the quenching of energy depositions in nuclear recoils, 
we apply correction factors obtained from previous stud¬ 
ies [3HH12]. Figure shows the energy spectrum of all 
crystals with energy depositions, after corrections factors 
are applied. 


B. Single-hit and multiple-hit events 

An important detail that is neglected in previous at¬ 
tempts UMllIS] to explain the DAMA signal modu¬ 
lation with muon-induced backgrounds is the acceptance 
for single-hit events. Figure shows the number of crys¬ 
tals in events with Apep > 2 keV per crystal, and in 
which at least one crystal in the event has a total en¬ 
ergy absorption of 2-6 keV. It is shown that < 9% of 


C. Events in the signal region 

In this section, we present the number of muon-induced 
single-hit events predicted by our simulation. The distri¬ 
bution of the energy deposited in crystals in single-hit 
events is shown in Figure For Apep < 20 keV, the 
muon-induced background is dominated by isolated neu¬ 
trons. 

In the range 2-6 keV there are 245 muon-induced 
events predicted over a period equivalent to twenty years. 
The total sensitive mass of the DAMA/LIBRA detec¬ 
tor is 242.5 kg, therefore we predict the rate of muon- 
induced events in this energy range to be 3.49 x I0“® 
counts / day / kg / keV with approximately 6% statisti¬ 
cal uncertainty. We estimate the systematic uncertainty 
to be approximately 30% by comparing different predic¬ 
tions of $((, as presented in Section [HT] 

We are able to compare our prediction to the conserva¬ 
tive estimate presented in Ref. HZlllH] which is in agree¬ 
ment with our results. 
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Energy deposition per crystal [keV] 


FIG. 4. The distribution of the total energy deposition in 
single-hit events. The blue dashed line and the black solid line 
show the sum of all energy depositions before and after energy 
resolution is considered, respectively. The stacked colored 
bars indicate the relative fraction of all events (before energy 
resolution is considered) attributed to events in which only 
neutrons deposit energy (yellow) and other events (blue). The 
equivalent of twenty years of muon-induced data is presented. 

The calculated event rate accounts for ~ 0.3% of the 
modulation amplitude reported by DAMA of (1.12 ± 
0.12) X 10“^ counts / day / kg / keV [1]. It is clear 
from this comparison that, even if the systematic un¬ 
certainty is bigger than our estimates, no muon-induced 
background can be used to explain the observed signal 
modulation. Our simulations (Figures and also 
show that, if muon-induced backgrounds could explain 
the DAMA data, one should expect a non-negligible mod¬ 
ulation of the muon-induced background above 6 keV, as 
well as for events with multiple hits, which is not seen by 
DAMA. 


V. DISCUSSION 

In this section we will argue that muon-induced neu¬ 
trons cannot explain the DAMA data, even before any 
estimate of is performed. 

We start the discussion in a general way, by considering 
any possible source of modulated signal, including dark 
matter, as has been done in Ref. |13]. The measured rate 
of events at DAM A/LIBRA is clearly dominated by ra¬ 
dioactive background above 6 keV, which imposes a strict 
limit on any interpretation of the modulated signal. This 
radioactive background is almost flat at low energies [33] , 
with the exception of a peak from at about 3 keV, 
which agrees with the DAMA measurements. To pre¬ 
serve the shape (‘flatness’) of the radioactive background 
in the region 2-6 keV, the total signal should be small 
and hence, the modulated fraction of the signal should 
be large. As an example, the measured modulated signal 
rate of 0.019 counts / day / kg / keV at 2-3 keV, assumed 
to be 5% of the total (average) signal, will give the total 
signal rate of 0.38 counts / day / kg / keV. This is already 
a significant fraction of the total measured rate at 2-3 keV 
(about 30%), requiring the radioactive background rate 


to drop by 30% at this energy whilst maintaining a flat 
background above 6 keV. No model of radioactivity pre¬ 
dicts a dip in the background below 6 keV [13]. 

Let us now consider muon-induced backgrounds within 
this context. We assume that and 4)'^ are modu¬ 
lated in a similar way, linked to the mean muon en¬ 
ergy at LNGS [H]. The LVD and Borexino [33] 
experiments have observed a muon flux modulation in 
the range of 1.3-1.5% of the total If the modu¬ 
lated signal in DAMA is due to a muon-induced effect, 
then the total rate of this ‘effect’ will be 0.0112/0.014 « 
0.8 counts / day / kg / keV. This is approximately 
equal to the total rate of ~ 1 counts / day / kg / keV 
observed by DAMA in the 2-6 keV energy range [4]. 
The effect is more dramatic in the 2-3 keV energy 
range, where the modulated signal is approximately 
0.0190 counts / day / kg / keV [4]. This would im¬ 
ply a total muon-induced background of 0.0190/0.014 « 
1.4 counts / day / kg / keV, which is higher that the 
total rate of events observed by DAMA. This is excluded 
by radioactivity models [43] . 

It is clear from the latter discussion that for any ex¬ 
planation of the DAMA signal to be consistent with the 
measured spectrum of events, it must satisfy the follow¬ 
ing qualitative criteria: 

• The amplitude of the effect must be very small com¬ 
pared to the DAMA event rate. 

• The modulation amplitude of the effect must not 
be much smaller than the average amplitude of the 
effect. 

• Any effect not satisfying the latter two criteria im¬ 
plies that there is a new model of suppressed ra¬ 
dioactivity in the region 2-6 keV, that does not ap¬ 
ply above 6 keV. 

• The modulation of the effect must only affect 
single-hit events, whilst disregarding multiple-hit 
events. 

• The explanation must simultaneously predict the 
phase and the period of the modulation. 

An explanation which incorporates muon-induced back¬ 
grounds cannot satisfy these criteria. 

VI. CONCLUSIONS 

We have presented an accurate simulation of the muon- 
induced background in the DAMA/LIBRA experiment, 
in response to proposals to explain the observed DAMA 
signal modulation with muon-induced neutrons. We have 
performed a full simulation of the DAMA/LIBRA appa¬ 
ratus, shielding and detector housing using Geant4.9.6. 

We have calculated the muon-induced neutron flux in 
LNGS to be 4>0 = 1.0 x 10“® cm“^ s“^ (without back- 
scattering), which is consistent with previous Simula- 









5 


tions. After selecting events which satisfy the DAMA sig¬ 
nal region criteria, our simulation predicts a background 
rate of 3.49 x 10“® counts / day / kg / keV. This accounts 
for approximately 0.3% of the modulation amplitude. We 
find that one would expect a non-negligible modulation 
of muon-induced background above 6 keV, as well as for 


events with multiple hits, which is not seen by DAMA. 

We conclude from our study that muon-induced neu¬ 
trons are unable explain the DAMA data. Furthermore, 
a large signal event rate, independently of the source of 
this signal, is inconsistent with radioactive background 
models. 


[1] R. Bernabei et al.. The DAMA/LIBRA apparatus, 
Nuclear Instruments and Methods in Physics Research 
Section A: Accelerators, Spectrometers, Detectors and 
Associated Equipment 592 (2008), 297 - 315, 

arXiv:0804.2738 [astro-ph] 

[2] R. Bernabei et ah. Particle Dark Matter and 
DAMA/LIBRA, AIP Conf.Proc. 1223 (2010), 50-59 
arXiv:0912.0660 [astro-ph.GA] 

[3] R. Bernabei et ah. Dark Matter Investigation By 
DAMA At Gran Sasso, International Journal of Modern 
Physics A 28 (2013), 1330022, arXiv: 1306.1411 
[astro-ph.GA] 

[4] R. Bernabei et ah. Final model independent result of 
DAMA/LIBRA-phasel, Eur.Phys.J. C73 (2013), 2648, 
arXiv:1308.5109 [astro-ph.GA] 

[5] CDEX Collaboration, Limits on light weakly interacting 
massive particles from the CDEX-1 experiment with a 
p-type point-contact germanium detector at the China 
Imping Underground Laboratory, Phys. Rev. D 90 
(2014), 091701, arXiv:1404.4946 [astro-ex] 

[6] Z. Ahmed et al., Results from a Low-Energy Analysis of 
the CDMS II Germanium Data, Phys. Rev. Lett. 106 
(2011), 131302, arXiv:1011.2482 [astro-ph.CD] 

[7] E. Armengaud et ah. Search for low-mass WIMPs with 
EDELWEISS-II heat-and-ionization detectors, Phys. 
Rev. D 86 (2012), 051701, arXiv: 1207.1815 
[astro-ph.CO] 

[8] LUX Collaboration, First Results from the L UX Dark 
Matter Experiment at the Sanford Underground 
Research Facility, Phys. Rev. Lett. 112 (2014), 091303, 
arXiv:1310.8214 [astro-ph.CO] 

[9] SuperCDMS Collaboration, Search for Low-Mass 
Weakly Interacting Massive Particles with SuperCDMS, 
Phys. Rev. Lett. 112 (2014), 241302, arXiv: 1402.7137 
[astro-ex] 

[10] J. Angle et ah. Search for Light Dark Matter in 
XENONIO Data, Phys. Rev. Lett. 107 (2011), 051301 
arXiv:1104.3088 [astro-ph.CO] 

[11] E. Aprile et ah. Dark Matter Results from 225 Live 
Days of XENONIOO Data, Phys. Rev. Lett. 109 (2012), 
181301, arXiv:1207.5988 [astro-ph.CD] 

[12] K. Blum, DAMA vs. the annually modulated muon 
background, arXiv: 1110.0857 [astro-ph.HE] 

[13] J. Ralston, One Model Explains DAMA/LIBRA, 
CoGENT, CDMS, and XENON, arXiv: 1006.5255 
[astro-ph] 

[14] D. Nygren, A testable conventional hypothesis for the 
DAMA-LIBRA annual modulation, arXiv: 1102.0815 
[astro-ph] 

[15] S. Chang et ah. Statistical tests of noise and harmony 
in dark matter modulation signals, Phys. Rev. D 85 
(2012), 063505, arXiv: 1111.4222 [hep-ph] 

[16] E. Fernandez-Martinez and R. Mahbubani, The Gran 


Sasso muon puzzle, J. Cosmol. Astropart. Phys. 07 
(2012), 029. 5 p, arXiv: 1204.5180 [astro-ph.HE] 

[17] R. Bernabei et ah. No role for muons in the DAMA 
annual modulation results. The European Physical 
Journal C 72 (2012), arXiv: 1202.4179 [astro-ph.GA] 

[18] R. Bernabei et ah. No role for neutrons, muons and 
solar neutrinos in the DAMA annual modulation 
results. The European Physical Journal C 74 (2014), 
arXiv:1409.3516 [hep-ph] 

[19] J. H. Davis, Fitting the Annual Modulation in DAMA 
with Neutrons from Muons and Neutrinos, Phys. Rev. 
Lett. 113 (2014), 081302, arXiv: 1407.1052 [hep-ph], 

[20] H. Wulandari et ah. Neutron background studies for the 
CRESST dark matter experiment, arXiv:0401032 
[hep-ex] 

[21] R. Persian!, Measurement of the muon-induced neutron 
flux at LNGS with the LVD experiment, PhD. Thesis, 
Universita di Bologna (2011) 

[22] M. Selvi, Analysis of the seasonal modulation of the 
cosmic muon flux in the LVD detector during 
2001-2008., Proceedings of the 3P* ICRC (2009). 

[23] Borexino Collaboration, Cosmic-muon flux and annual 
modulation in Borexino at 3800 m water-equivalent 
depth, JCAP 1205 (2012), 015, arXiv: 1202.6403 
[hep-ex] 

[24] D. Barker, D.-M. Mei, and C. Zhang, Muon-induced 
background study for an argon-based long baseline 
neutrino experiment, Phys. Rev. D 86 (2012), 054001 
arXiv:1202.5000 [physics.ins-det] 

[25] S. Barbeau, P. et ah. Comment on “Fitting the Annual 
Modulation in DAMA with Neutrons from Muons and 
Neutrinos”, Phys. Rev. Lett. 113 (2014), 229001, 
arXiv:1409.3185 [hep-ph] 

[26] S. Agostinelli et ah, Geant4 - a simulation toolkit, 

Nucl.Instrum.Meth 506 (2003), 250 - 303 

[27] V. Kudryavtsev, Muon simulation codes MUSIC and 
MUSUN for underground physics, Computer Physics 
Communications 180 (2009), 339 - 346, 

arXiv:0810.4635 [physics.comp-ph] 

[28] P. Antonioli, C. Ghetti, E. Korolkova, V. Kudryavtsev, 
and G. Sartorelli, A three-dimensional code for muon 
propagation through the rock: MUSIC, Astroparticle 
Physics 7 (1997), 357 - 368, arXiv: 9705408 [hep-ph] 

[29] We define ‘primary’ particles as those which are present 
at the surface at the Earth, and ‘secondary’ particles as 
those which are subsequently produced. 

[30] Cross Section Evaluation Working Group, ENDF/B-VII 
data files may be accessed through the national nuclear 
data center, www.nndc.bnl.gov . 

[31] H. Araujo et ah, Measurements of neutrons produced by 
high-energy muons at the Boulby Underground 
Laboratory, Astroparticle Physics 29 (2008), 471 - 481, 
arXiv:0805.3110 [hep-ex] 




6 


[32] Borexino Collaboration, Muon and cosmogenic neutron 
detection in Borexino, Journal of Instrumentation 6 
(2011), P05005, arXiv: 1101.3101 [physics . ins-det] 

[33] B. Schmidt et ah. Muon-induced background in the 
{EDELWEISS} dark matter search, Astroparticle 
Physics 44 (2013), 28 - 39, arXiv: 1302.7112 
[astro-ph.CO] 

[34] L. Reichhart et ah. Measurement and simulation of the 
muon-induced neutron yield in lead, Astroparticle 
Physics 47 (2013), 67 - 76, arXiv: 1302.4275 
[physics.ins-det] 

[35] H. Wulandari et ah. Neutron flux at the Gran Sasso 
underground laboratory revisited, Astroparticle Physics 
22 (2004), 313 - 322, arXiv:0312050 [hep-ex] 

[36] We define the back-scattered neutron flux as including 
independent counts from neutrons that re-enter the 
cavern due to scattering in the surrounding rock. 

[37] A. Ferrari et al. FLUKA: A multi-particle transport 
code. CERN, Geneva. arXiv: 1406.6081 
[astro-ph.IM] and references therein. 2005. 

[38] D. Tovey et al.. Measurement of scintillation efficiencies 
and pulse-shapes for nuclear recoils in Nal(Tl) and 
CaF-2(Eu) at low-energies for dark matter experiments, 
Phys.Lett. B433 (1998), 150-155 


[39] N. Spooner et ah. The Scintillation efficiency of sodium 
and iodine recoils in a Nal(Tl) detector for dark matter 
searches, Phys.Lett. B321 (1994), 156-160 

[40] H. Chagani et ah. Measurement of the guenching factor 
of Na recoils in Nal(Tl), Journal of Instrumentation 3 
(2008), P06003, arXiv:0806.1916 [physics . ins-det] 

[41] G. Gerbier et ah. Pulse shape discrimination with 
Nal(Tl) and results from a WIMP search at the 
Laboratoire Souterrain de Modane, Astropart.Phys. 11 
(1999), 287-302 

[42] E. Simon et ah, SICANE: A Detector array for the 
measurement of nuclear recoil guenching factors using a 
monoenergetic neutron beam, Nucl.Instrum.Meth. A507 
(2003), 643-656, arXiv:0212491 [astro-ph] 

]43] V. Kudryavtsev, M. Robinson, and N. Spooner, The 
expected background spectrum in Nal dark matter 
detectors and the DAMA result, Astroparticle Physics 
33 (2010), 91 - 96, arXiv: 0912.2983 [hep-ex] 

[44] V. Kudryavtsev, N. Spooner, and J. McMillan, 
Simulations of muon-induced neutron flux at large 
depths underground, Nuclear Instruments and Methods 
in Physics Research Section A: Accelerators, 
Spectrometers, Detectors and Associated Equipment 
505 (2003), 688 - 698, arXiv:0303007 [hep-ex] 


